Reply to the Reply to the Comment on "Z 2 -slave-spin theory for strongly correlated 

fermions" 



Alvaro Ferraz 1 , Evgeny Kochetov 1,2 
1 International Institute of Physics - UFRN, Department of Experimental and Theoretical Physics - UFRN, Natal, Brazil and 
2 Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 14-1980 Dubna, Russia 



The authors of the Reply [l| to our Comment claim 
that the equation 

Ai\phys) = i[A h H]\phys) = DiA^phys) = (1) 

insures that the on-site operator constraint A = se- 
lects the physical subspace consistently with the time 
evolution generated by the Hamiltonian H. Here D t is 
some nonlocal operator constructed out of the spin and 
fermion operators. This is indeed the case since an ex- 
plicit solution of Eq.fl]) reads A(^) = exp (IDj)Aj. As 
AAphys) = 0, one gets Ai(t)\phys) — at any given time 
instant. 

No calculations are in fact needed to see this, since the 
equation Ai(t)\phys) = e~ ltH A i e' ,tH \phys) = follows 



mutes with H n 



the constraint Qi = must now be 
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simply from the observation that the Hamiltonian H does 
not couple the physical and unphysical subspaces. This 
occurs because of the fact that it commutes with the pro- 
jection operator Qi := A? which has eigenvalues and 1 
corresponding to the physical and unphysical subspaces, 
respectively. In case the dynamics never takes the state 
outside the physical subspace, the requirements Aj = 
and Qi = are truly equivalent. 

However, as soon as any approximate treatment that 
couples the physical and unphysical spaces is applied, 
those requirements are no longer identical. In this case 
Qi = remains the only constraint available to select the 
physical subspace. Since the operator Qi no longer corn- 



enforced explicitly at any given time instant. This can 
be done by adding to the Hamiltonian a global Lagrange 
multiplier in the form AJ^Qi. Since the eigenvalues 
of the on-site projection operator Qi are either or 1, 
the limit A — > +oo singles out the physical subspace. Q 
Only after this is done, can one safely use the equation 
Ai = 0. Within a mean-field (MF) approximation, one 
may however replace the operator equation Qi = with 
its ground-state expectation value, (Qi) = 0. In spite of 
that it still eliminates the unphysical states. In contrast, 
the constraint expectation value (A) can be nullified by 
the unphysical subspace as well. [2j] Consequently, a MF 
theory based solely on the condition (Ai) = becomes 
inconsistent since uncontrolled unphysical degrees of free- 
dom contribute to the theory as well. 

To conclude, Reply [lj does not refute the criticism 
exposed in our Comment 0|. 
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[3] This trick cannot be used to impose the requirement At = 
0, since the operator Ai has an eigenvalue —1. 
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The authors of ref. [l[ introduce a new slave-spin (SS) 
mean-field (MF) approach to the Hubbard model of 
strongly correlated electrons in terms of Ising spins and 
standard fermion operators, by constructing an enlarged 
Hilbert space with both physical and unphysical states. 
The purpose of this construction is clear: this formulation 
involves a much smaller number of slave fields, compared 
with more standard slave particle frameworks. In this 
way the complexity of the original model is reduced to a 
more readily tractable system of Ising spins in a trans- 
verse field coupled to free fermions. In this Comment, we 
show that the unphysical states are not properly excluded 
and despite the claimed good agreement with previous 
results, the consistency of their MF theory is unjustified. 

The paper [l[ may be considered as an extension and 
an improvement of the formalism developed earlier in [2j . 
It can also be thought of as a minimal formulation of 
previous slave-spin representations. Q 

The single-band Hubbard model is written in the form 



U 



(1) 



The hopping amplitude is denoted by tij and U is the 
onsite repulsion. The operator c^J creates an electron at 
site i with spin a, and rij = c\ Ci a . 

The authors introduce auxiliary local pseudospin vari- 
ables Ii = (If , If , If) with eigenstates If|±)i = ±||±)j. 
The physical electron operator is then represented as 



it) 



O TX j(t) 

A1 i Jia i 



(2) 



where is an auxiliary fermion operator defined on the 
space spanned by the vectors |0)j, | f)i, I I t4-)i- [H 

The SS representation ([2]) expands the Hilbert space. 
The authors identify the on-site physical Hilbert space to 
be spanned by the states 

n phys = {\+)\o), |-)|t), |-)ll), l+)IU)}- (3) 

The unphysical states are then 

n unphys = {|->|o>, |+)| t>, l+)l I), |-)l tl». (4) 

To obtain a faithful representation of the original prob- 
lem, the authors claim to single out the physical subspace 
by the local constraint A := If + 1/2 — (n, — l) 2 = 0. 
This equation justifies the substitution in (1) 



I) 2 



1/2), 



(5) 



which appears as a necessary step to recover the trans- 
verse field Ising model representation. 
The projection of the Hubbard model 



H' 



ija 
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1. 



onto the physical subspace is equivalent to the original 
model (fTJ). Such a projection can be achieved explicitly 
by imposing the requirement Qi :— = 1/2 + If[l — 
2{rii — l) 2 ] = 0, for each lattice site. The projection op- 
erator Qi commutes with Hamiltonian (flj and generates 
a U{1) local gauge symmetry which takes care of the 
redundancy of the decomposition @. However, notice 
that the operator constraints A; = and Qi = are not 
equivalent to each other. 

The, so far, exact SS formulation is then treated at MF 
level by fully decoupling the SS and the auxiliary fermion 
operators in the form H' — > H MF — Hf + Hj, where 



— 9ij tij fia fja , 



(7) 



Here the hopping amplitude and the Ising exchange cou- 
pling are renormalized by the factors gij and \iji re_ 
spectively. These are to be determined through the self- 
consistency equations 

9ij = i(IfIf) Hl , Xij = 4E((4/ j .)h / + cc). (8) 



To proceed, the authors replace the constraint Aj = 
by its ground state expectation value, ignoring com- 
pletely the condition Qi = 0. In the enlarged 
Hilbert space, the Qi operator is a projection matrix, 

Qi\T~Lphys)i — 0, Qi\Wunphys)i — \'Hunphy s) i ■ Because of 

this, the local MF constraint (Qi) = still eliminates the 
unphysical states. In contrast, the same condition does 
not hold for the constraint expectation value (Aj) that 
can be nullified by the unphysical subspace as well. To 
see this, let us pick up a state 



i$) = cf ; i-) i |o) i +cn-) i |n) 



(2) 



+cf ) \ + ) l \ t)< + cH+)i| i)i G Uunphys 



.(4) 



2 



where |c| 1} | 2 + \cf ] \ 2 = {cf^ 2 + | Cl (4) | 2 . It then follows 
that = 0. Consequently, the condition (A t ) = 

no longer discriminates between physical and unphysical 
states (as opposed to the requirement (Qi) = 0), and 
the resulting SS MF theory becomes inconsistent. The 
substitution (5) now implies that the uncontrolled un- 
physical degrees of freedom contribute to the theory as 
well. 

The authors refer to their MF approach as a Z2 gauge 
theory. They claim that the local Z2 transformations, 
fi -> UifiUi = -fi, If -t UilfUi = -If, Ui = 1 - 2Qi, 
respect the MF decomposition and "in this sense, the 
mean- field ansatz breaks the U(l) gauge theory down to 
Z2" -[ll This statement may cause some confusion. [8| It 
is of course legitimate to restrict the space of variational 
states to a subspace generated by the states related by 
the Ui transformations. However, this does not result 
in a Z2 invariant MF theory: the gauge-related varia- 
tional states still belong to the different Hamiltonians, 
H MF and UiH MF Ui ^ H MF . In fact, a Z2 gauge invari- 
ant MF theory only arises if one considers both functions 
Xij and gij as Hubbard-Stratanovich dynamical fields, 
Xij = XijVij, 9ij = SijVij- Here Xij and are precisely 
the solutions to the saddle point Eqs.([8]), and cry = ±1 
is an Ising dynamical gauge field that takes care of the 
Z2 gauge fluctuations beyond the saddle point approxi- 
mation. 
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